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' Abstract 
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Efficient algorithms are presented for constructing spanners in geometric intersection graphs. 
For a unit ball graph in M k , a (1 + e)-spanner with 0(ne~ k+1 ) edges is obtained using efficient 
partitioning of the space into hypercubes and solving bichromatic closest pair problems. The 
spanner construction has almost equivalent complexity to the construction of Euclidean min- 
C > \ imum spanning trees. The results are extended to arbitrary ball graphs with a sub-quadratic 

running time. 

For unit ball graphs, the spanners have a small separator decomposition which can be used 
to obtain efficient algorithms for approximating proximity problems like diameter and distance 
queries. The results on compressed quadtrees, geometric graph separators, and diameter 



o 
u 

approximation might be of independent interest. 

1 Introduction 

> 

; Let G = (V, E) be a weighted graph, and let de(u,v) be the length of shortest path between 

vertices u and v in G. For any fixed e > 0, a stretch (1 + e)-spanner of G is a subgraph G' 
such that for all pairs of vertices u and v, dc' (u, v)/dG(u, v) < (1 + e). Spanners are important 
structures because they provide a way for approximating the graph in an economical way. Spanner 
constructions have been widely investigated for general graphs, and complete Euclidean graphs, 
also with additional properties like weight, diameter, degree ACIM99] ICK951 IADM+95j . 

We present a new method for producing spanners of geometric graphs based on a hierarchical 
decomposition of the plane into tiles of various sizes. Our constructions are also more general, 
^ as they are not restricted to complete Euclidean graphs, but extend to geometric unit disk and 

arbitrary disk graphs, as well as their higher dimensional versions, the ball graphs. In all cases, 
edge lengths are given by Euclidean distances, but not all edges have to be present in our graphs. 
The difficulty in constructing a spanner for the disk graph metric when compared to the metric 
induced by a complete Euclidean graph is that two points that are close in space are not necessarily 
close under the graph metric. 

Intersection graphs are graphs whose vertices are represented by sets such that two vertices are 
adjacent if and only if the corresponding sets have a non-empty intersection. A disk graph is an 
intersection graph of disks in the plane. We consider weighted disk graphs where the weight of an 
edge is the Euclidean distance between centers. Such graphs have been used widely to model the 
communication between objects in VLSI |MC80| and recently in the context of wireless networks 
KMR0ll lLi03| . In the wireless networks they model the fact that two wireless nodes can directly 
communicate with each other only if they are within certain distance. Many restricted versions 
of disk graphs like unit disk graphs, quasi unit disk graphs, A-precision unit disk graphs have also 
been investigated as models for wireless networks |Li03l IKZ031 lHMR + 9"8 . 
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Spanners are important for disk graphs because restricting the size of a network reduces the 
amount of routing information. One of the requirements of topology control is to produce a low 
stretch spanner |Raj02[ fLi03j . Li et al. LWW01 showed that deleting certain edges from the 
Yao-graph (introduced in |Yao82j ) results in a bounded stretch and bounded degree spanner for 
the special case where all disks have unit radius. For solving proximity problems Gao and Zhang 
GZ05 gave a construction of a c- well-separated pair decomposition (introduced in |CK95j ) for 
unit disk graphs with O(nlogn) pairs, for any constant c > 1 . They also show that in M k at least 
0(u 2 ~ 2 / k ) pairs are needed. However, for general disk graphs no sub-quadratic well-separated 
pair-decomposition is possible. 

2 Preliminaries 

Let V be a set of points in M k for any fixed dimension k. Let T> be a set of n balls such that (i) D p 
is centered at p £ V, and (ii) D p has radius of r p . Balls D p and D q intersect if d(p, q) < (r p +r q ), 
where d(., .) denotes the Euclidean metric. The disk graph G is a weighted graph where an edge 
e = (p, q) with weight d(p, q) exists if D p and D q intersect. Let dG denote the metric induced by 
the connected graph G on its nodes by shortest paths. We re-scale the balls such that the largest 
radius equals one. The global stretch factor of T> is then defined as p{T>) = 1 /min{r p : D p £ T>}. 

In the plane our algorithms use a variant of quadtrees. For a node t, denote by P(t) the parent 
of t in the tree. A point (x,y) is contained in a node t representing a square with center (x t) yt) 
and length l t in the quadtree iff x t — lt/2 < x < x t + l-t/2 and y t — l-t/2 < y < Vt + l-t/2. For a 
set of squares T in the quadtree a point is contained in T iff there exists t S T, such that point is 
contained in t. 

To avoid ambiguities, throughout the paper we refer to the vertices of a graph as vertices and 
vertices of a tree as nodes. We assume w.l.o.g. that e _1 is a power of 2. Floors and ceilings are 
omitted throughout the paper, unless needed. Note that starting with a fixed e we get a (1 +ce) 
spanner for a fixed constant c. For simplicity we describe all the algorithms for k = 2 and then 
state the generalizations to higher k. We use the notation 0(f) = 0(f poly log f). 

2.1 Modified Yao-graph and Our Results 

We use a modified version of the Yao-graph construction. Let C(p) = {c-\ (p), . . . , c e -i (p)} be a 
collection of e _1 cones such that (i) each cone has apex at p E V, (ii) each cone has an opening 
angle of 2ne, and (iii) the union of these cones covers R 2 . We define a modified Yao-graph Y* in 
the following way. The vertices of Y are the points of V . For each p 6 V and 1 < i < e _1 add 
a directed edge from p to point q contained in Ct(p) if: (i) q is closest to p in c-t(p), (ii) r q > r p , 
and (iii) the edge (p , q ) exists in G . We use Y to denote the undirected counterpart of Y . 

Lemma 1. Let (u,v) he an edge in the disk graph G. Then there exists a path in the spanner Y 
such that dv(u,v) < (1 + e)d(u,v). 

Proof. Let r u < r v . Consider the cone Ct(u) containing v. Let (u,w) be the edge added in 
Y from u to w contained in Ci(u). The distance d(w,v) < d(u,v) — (1 — 0(e))d(u,w). Since 
r w > i"ui we know that there exists an edge (w,v) in G. The proof can be completed by using 
induction over distances. □ 

From the above lemma by summing over all edges of a path in G we also get that Y is a (1 + e) 
spanner for the disk graph G. For ball graphs in M k , Y has 0(ne~ k+1 ) edges |Luk99j . For unit 
ball graphs the modified Yao-graph is same as the original Yao-graph with long edges removed. 
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Figure 1: The quad-dissection procedure for a set of points and the corresponding forest V. The 
donut shaped nodes are the interesting nodes in V. 

In the remaining paper we improve upon this result in the following manner. All constructions 
are for (1 + e)-spanner. 

Unit Ball Graphs: For unit ball graphs in M k we solve bichromatic closest pair problems to pro- 
duce spanners with 0(ne _k+1 ) edges. A spanner can be constructed in time 0(nlogne _1 +ue~ 2 ) 
for k = 2 by solving post-office problems. In higher dimensions using the currently best algorithm 
of Aggarwal et al. AESW91] for solving the bichromatic closest pair problem, we can construct 
a spanner in expected running time 6(n 4 / 3 e~ 2 + ue~ 3 ) for k = 3 and 0(u 2 ~ 2 /T k /21+i)+6 e -k+l _|_ 
ne _k ) for k > 4, where 6 is any positive constant. In addition, we show that constructing (1 + e)- 
spanners for unit ball graphs has almost equivalent randomized complexity to solving bichromatic 
closest pair problem for n points. Additionally the spanners for unit ball graphs have small ge- 
ometric separator (hereditary), which we exploit to obtain faster approximation algorithms for 
proximity problems like diameter and distance queries. 

Ball Graphs: We generalize the above result to ball graphs in M k . In this case we replace bichro- 
matic closest pair problems by halfspace range searching problems. The spanner has 0(ue~ k+1 ) 
edges. For the interesting case when p{T>) is polynomially bounded by using the currently best 
algorithm of Aggarwal and Matousek |AM95j for half space range searching, we obtain a running 
time 0(n 2d+6 e~ kd + ne~ k ), where d = 1 — 1 /( [k/2J + 2) and 6 is any positive constant. 

3 Spanners for Unit Disk Graphs 

To illustrate our key ideas we first describe an algorithm that produces a spanner G ' of unit disk 
graph G with 0(ne~ 2 ) edges. Later we describe the procedure for sparsifying the spanner. We 
first describe a high level idea of our algorithm and prove the claimed stretch factor, and then 
define an algorithmic version of the construction. The spanner G ' is constructed in the following 
manner. 

Quad-Dissection: The input to our algorithms is a set of points (V) in M 2 . The bounding box of 
V is a smallest rectangle enclosing the points and we assume the left bottom corner as the origin. 
An L-grid is defined by horizontal and vertical line segments drawn at y = LZ and x = \JL within 
the bounding box. Our spanner construction involves recursively partitioning the plane using 
a simple variant of quadtrees. A quad-dissection of the the L-grid is a recursive partition into 
smaller squares. We view it as a 4-ary forest with the root nodes as the non-empty squares in the 
L-grid. Each square is partitioned into four equal squares, which form its children. We continue 
partitioning the non-empty squares until each disk center is contained in a separate square. 
Constructing the Forest: Let V denote the forest from the quad-dissection of the e-grid. Roots 
represent the non-empty squares of the e-grid. V is a collection of disjoint trees, each of which 
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Close Neighborhood: N c (t) 
Far Neighborhood: N (t) 



Figure 2: Close and Far neighborhoods of a node t in Roots. Concentric circles are 0(e) apart. 

is rooted at a node of Roots. Note that the set of nodes at depth I in the forest corresponds to 
the set of non-empty squares defined by the e2~ l -grid. Define C(t) as set of disk centers lying 
contained by the node t. A node t of the forest is called interesting if the set of disk centers 
contained it is not equal to the set of disk centers in its parent, i.e., C(t) 7^ C(P(t)). By definition 
all nodes in Roots are interesting. It immediately follows that V has at most 2n interesting nodes. 
Choosing the Representatives: For every leaf node t G V we chose the disk center in C(t) as 
its representative R t . For every internal node t G V we pick one of its children t' and set R t to 
R t ,. " 

Neighborhood of Nodes: For every t G Roots, define its close neighborhood Nc(t), as the set 
of all nodes in Roots which are within a distance of 2 — \/2e. The distance between squares can 
be thought as distance between their centers. For every t G Roots, define its far neighborhood 
Np(t), as the set of all t' G Roots which are within distance (2 — \/2e,2 + V2~e] and for which 
there exists at least one edge between vertices of C(t) and C(t') in G. For every interesting node 
t ^ Roots, define its close neighborhood Nc(t), as the set of all nodes at depth dt which are 
within a distance of 2~ dt from t. For every interesting node t ^ Roots set its far neighborhood 
N F (t) to 0. 

Edges in Spanner G': For every interesting node t, we add a edge between R t and R t ', where 
t' G Nc(t). Additionally for every node t G Roots, we put an edge in G between disks C(t) and 
C(t') into the spanner, where t' G Np(t). 

Lemma 2. The spanner G' has 0(ue~ 2 ) edges and is a subgraph of the unit disk graph G. 

Proof. The nodes at depth dt in V are defined by the set of non-empty squares of e2~ dt -grid. 
Thus the set Nc(t) is of size 0(e~ 2 ) and for every interesting node we add 0(e~ 2 ) edges from its 
representative R t . There is an edge (Rt, Rt') with t' G Nc(t) in G because the Euclidean distance 
between them is at most 2. 

Also for every t G Roots, the nodes in Np(t) are the nodes that lie between two concentric 
circles 2s/2e apart. Since nodes in Roots are defined by the e-grid, there are only 0(e _1 ) nodes 
in Np(t). Thus we only add 0(e _1 ) additional edges for every node in Roots and these edges are 
added only if they are present in G. Since there are at most 0(n) interesting nodes, we get the 
claimed result. □ 



Lemma 3. Let u and v be disk centers, such that there exists a node t in V, with u,v G C(t). 
Then there exists a path in the spanner G' such that dc' (u, v) < (1 + c u e)d(u, v) for some constant 
c u . 

Proof. The maximum square length is e. This implies that if there exists a node t with 
u,v G C(t), then dc(u,v) = d(u,v). The proof is by induction over the size of C(t). The base 
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case of |C(t)| = 1 is trivial. Now in the subtree rooted at t, consider the maximum depth such 
that there exist nodes a and b with u G C(a), v G C(b) with an edge between R Q and Rb in G'. 
There exists such a and b because interesting nodes with the same parent have an edge between 
their representatives. The case where u = R a or v = Rb is simpler and for now we assume u/R a 
and v / Rb- 

In the subtree rooted at a, find the closest descendant a' of a which is interesting with 
u G C(a')- Let b' be the node at same depth as a' with v G C(b'). There exists no edge between 
between representatives of a' and b' (because a and b are the deepest nodes with this property). 
This implies that d(u,v) > 2~ da '. Also if not a', at least P(a') contains both u and R a . Therefore 
for some constant c s , we get d(u, R a ) < c s e2~ d a' and thus d(u, R Q ) < c s ed(u, v). Using similar 
arguments one can show that d(Rb,v) < c s ed(u,v). 

By the inductive hypothesis we know, dG'(u-,R a ) < C +c u e)d(u,R a ) and dG'(Rb> v ) < (1 + 
c u e)d(Rb, v). In G' there exists a path from u to v of length (1 + Cu.e)d(u, R Q ) + d(R a , Rb) + (1 + 
c u e)d(Rb,v). Finally putting together everything we get 

d G '(u ( v) = (1 +c li e)d(u,R Q ) + d(R a ,Rb) + (l + c u e)d(R b) v) 

< (2 + c u e)d(u, R Q ) + d(u,v) + [2 + c u e)d(R b , v) < (1 + c u e)d(u,v) 

The second inequality comes from applying the triangle inequality and the final step involves 
substituting d(u, R a ) and d(Rb,v) in terms of d(u,v). The constants satisfy c s < 2\/2 and 
Cu>4c s /(1 -2c,e). □ 

Lemma 4. Let (u,v) be an edge in the unit disk graph G. Then there exists a path in the spanner 
G' such that d.Qi{u,v) < (1 + c u e)d(u,v) for some constant c u . 

Sketch of Proof. If there exists a node t such that u, v G C(t), then by Lemma El the claim is 
true. Otherwise, consider the deepest nodes a, b in F, such that (i) u G C(a), v G C(b), and (ii) 
edge between u' G C(a) and v' G C(b) is in G'. Such a and b exist at least as nodes in Roots. 
As in LemmaElone can argue that d(u,u') < c s ed(u,v) and d(v',v) < c s ed(u,v). From Lemma 
121 we also know that dG'(u,u') < (1 +c u e)d(u,u') and dQ'[v',v) < (1 + c u e)d(v',v). The proof 
can be completed by applying the triangle inequality. □ 

Sparsifying the Spanner: Let C[t) be the collection of e _1 cones such that (i) each cone has 
apex at Rt, (ii) each cone has an opening angle of 2ne, and (iii) union of these cones covers the 
entire area represented by Nc(t). 

Instead of adding edges to every R t ', t' G Nc(t), we include the shortest edge between R t 
and R t ' in each cone. Since there are only e~ 1 cones, we conclude that the number of edges in 
the spanner is 0(ne _1 ). Using the analysis of Lemma ^along with Lemma |1J it follows that the 
graph is a (1 + e)-spanner. Let G' denote this sparsified spanner. 

3.1 Algorithmic Version 

We now describe the algorithmic counterparts for each step in the construction of G'. 
Constructing the Forest: We construct a compressed forest (compressed quadtrees) V, in 
which we only introduce the interesting nodes of V and shortcut the degree one internal nodes. For 
the sake of construction we introduce some terminology. Given two binary numbers X = ^ *h2 x 
and Y = ^\)i2 l , define a bitwise shuffle operation of X and Y as: 

X|Y=V^(2x i + -y i )4 i 
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Define agree(X, Y) as the smallest I such that for X / Y we have [X2 2l+1 J / [ Y22l+1 J- Note that 
if X|Y< X'|Y', then either agree(X,X') < agree(Y,Y') and X < X' or agree(X,X') > agree(Y,Y') 
and Y< Y'. We also define agree(X, Y, I'), with X / Y as the smallest I satisfying |_X2 2l+l ' +1 J 7^ 
LY2 2l+l ' +1 J. Define rouud(X, I) = [X2 l \/2 l . 

For a point p, let x(p) and y(p) denote its x and y-coordinates and let p s = x(p)kj(p). Sort 
the points in V in increasing order of their shuffled representations. Let V s = {pf ,p|, • • • , PiU 
denote this sorted array. Construct a list A s whose i th entry at = agreefpf.p^plog e -1 ). Start 
the forest with a different leaf node for every disk. We maintain a stack with top representing the 
last inserted entry and a special bottom of stack entry (—00, Null). When we scan a new element 
at in the list A s we do: (i) if at > atop, push (at,tt) onto the top of the stack and increment i. 
(ii) else, pop (a t op,"ttop) and ( Qtop-1 , ttop-1 ) and create a new node t in V with its children as 
t t0 p and t top _-|. Label edges (t,t t0 p) and (t,t top _-|) with a t0 p- Push (a t0 p-i)"t) on the stack. 

The depth of a node t, (d t ) is the defined by the label of the edge connecting it to its parent. 
The V so generated is the binary tree representation of the compressed quadtree. It can easily be 
converted into a compressed quadtree by just shortcutting the nodes which have the same depth 
as their children. The bottom left hand corner of t has coordinates (roun,d(x(R t ), d t + log e -1 ), 
rouud(u(R t ), d t + log e -1 )). Roots is the set of nodes at depth 0. 

The time for constructing the forest is O(ulogu), assuming agree 1 operation can be imple- 
mented in constant time. The shuffle operation needn't be performed explicitly as agree can be 
used to generate V s - 

Finding the close neighborhood: Finding the close neighborhood for nodes in Roots is straight 
forward. Because we work with V we construct a close pseudo-neighborhood (N^(t)) for a node t. 
The close pseudo-neighborhood consists of a set of nodes, where each element r of the set either 
belongs to Nc(t) or in V the node r is the maximal depth ancestor of a node in Nc(t). 

Now assuming we have found N^(t), we describe a method of constructing N^(t') for a child 
t' of t in F'. We access all the nodes in N^(t). For every a G N^(t), we look at the subtree 
rooted at a. In this subtree we do a BFS traversal with a modification that the subtree for any 
b G is accessed only if db < d t ' and Rj, is at most distance (1 + \fle)2~ dh away from R t /. The 
nodes at which the traversal ends and whose representatives are at most (1 + V2e)2~ d t' away 
from R t / define N^(t'). 

Any node g G V accessed during the traversal has a node f G V such that: (i) f is an ancestor 
of t' in P, (ii) d g = df, and (iii) the distance between f and g is 0{2~ de ). We charge the access 
to node g to node f (f needn't be in V). For the entire algorithm we can charge all the accesses 
to node g to different nodes which are at the same depth as g and only 0(2~ dg ) away (note there 
are only 0(e ) such nodes). This implies that the time for finding close pseudo-neighborhoods 
for all nodes is 0(ne~ 2 ). Alternatively, one can add edges to only the closest representative in 
each cone within the same time bound. 

Finding the far neighborhood: The far neighborhood of a node t G Roots is computed by first 
determining the set of nodes F(t) C Roots which are at distance (2 — \/2e, 2 + \[2e\ away from t. 
For each t' G F(t), we solve the bichromatic closest pair problem with one set as C(t) and other 
as C(t'). If the disks corresponding to the closest pair intersect, we add the corresponding edge 
into the spanner. It is well known that bichromatic closest pair in M 2 can be solved in O(nlogn) 
by using post-office problems. Since each set C(t) takes part in only 0(e _1 ) bichromatic closest 
pair problems, we get a time of 0(ulogne _1 ) for finding all far neighborhoods. 
We can summarize these results as follows. 

lr The authors believe that the same running time can be obtained by just using floor of binary logarithm (a 
simple operation) . Also the floor operation used here is not dangerous as we don't multiply variables. 
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Theorem 1. There exists an algorithm for producing a (1 + e) -spanner for unit disk graphs. The 
spanner has 0(ne~ 1 ) edges and can be computed in time 0(nlogrie~ 1 + ue~ 2 ). 

3.2 Extension to Unit Ball graphs 

The quadtrees becomes octrees in M k . The close neighborhood of any node is of size 0(e~ k ), 
but again one can sparsify to get it to 0(e~ k+1 ). The far neighborhood of any node in Roots 
is 0(e~ k+1 ). Therefore the total number of edges is 0(ne~ k+1 ). The bottleneck in the running 
time comes from solving bichromatic closest pair problems for finding far neighborhoods of nodes 
in Roots. We now state a general result describing the upper bounds of the spanner construction 
in terms of the upper bounds for bichromatic closest pair problem. 

Let f(n) be any function satisfying logf(u) = o(logn), i.e., f grows slower than polynomial. 
We assume in the rest of discussion that n c f (n) is an upper bound of on the time for computing 
a bichromatic closest pair for a total on n points, for some constant c > 1 . 

Lemma 5. For unit ball graphs in M k , the worst case running time for finding the far neighborhood 
for all the nodes in Roots is 0(n c f (n)e~ k+1 ). 

Proof. For every node t in Roots, we solve bichromatic closest pair problems between C(t) and 
C(t'), t' £ Roots and |C(t')| < |C(t)|. Since t participates in 0(e~ k+1 ) bichromatic closest pair 
problems, the cost of solving these problems (charged to t) is 0(|C(t)| 2c f (|C(t)|)e~ k+1 ). Since 
every disk center is contained in exactly one node of the Roots, the total cost is 



0( Y. IC(t)| c f(|C(t)|)]e- k+1 ) < 0(n c f(u)e- 



-k+1 - 
te Roots 



□ 

The currently best algorithm of Aggarwal and Matousek [AESW91 for finding bichromatic 
closest pair between P red and Q blue points in M k runs in expected time 0((PQ log P log Q) 2 ^ 3 + 
P log 2 Q + Q log 2 P) for k = 3 and in expected time 0((PQ) 1_l /(rV2l+1)+s + p i og Q + Q log P) for 
k > 4, where 5 is any positive constant. We can summarize the results for spanners in unit ball 
graphs as 

Theorem 2. There exists an algorithm for producing a (1 + e)- spanner for unit ball graphs in M k 
(k > 3). The spanner has 0(ue~ k+1 ) edges and can be computed in time 0(n c f (n)e — k+1 +ue~ k ). 

Corollary 1. The spanner G' can be constructed in expected time 0(n 4 / 3 e~ 2 + ne~ 3 ) for k = 3 

and for k > 4 it can be constructed in expected time 0(u 2 ~ 2 /T k /21+l)+6 e -k+i +ne~ k ), where 8 is 
any positive constant. 

The above result shows that finding spanner of unit ball graphs in M k is not much harder than 
computing bichromatic closest pair for n points in M k . In the other direction from the results 
of Eppstein Epp04 and Chan jCha99j we know that the randomized expected time bounds for 



constructing even the spanning forest of unit balls in M k and a Euclidean minimum spanning tree 
(or a bichromatic closest pair) in M k are within constant factors. To make the relation between 
the running time for spanner construction and bichromatic closest pair problem more precise we 
introduce a few definitions. Define exponent of a problem A with respect to input size n to be 

inf{c | there exists an algorithm for solving A with running time 0(u c )} 

Since the bichromatic closest pair and problem and Euclidean minimum spanning tree have asymp- 
totically same complexities AESW91, KLN99], their exponents are also equal. 
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Theorem 3. The exponent of (1 + e) - spanner fom unit balls in IR k is same as the exponent of 
bichromatic closest pair for n points in M k . 

Proof. From Theorem [51 we know that exponent of the (1 + e)-spanner is not more than the 
exponent of the bichromatic closest pair. For the other direction, once we have the spanner any 
graph traversal can be used to construct the spanning forest of unit balls. Therefore, from results of 



Eppstein Epp04 about randomized equivalence between spanning forest and bichromatic closest 
pair, we also get that the exponent of bichromatic closest pair is not more than the exponent of 
(1 + e)-spanner. □ 

4 Separators in Spanner Graph 

In this section we show that the spanner G' has a small geometric vertex separator. Our results 
can be seen as analogous to results by Smith and Wormald SW98] for separators in spanners 
for complete Euclidean graphs. We show that G' has an n 1— 1//,c e — k+1 / 2 - vertex separator, whose 
removal leaves two sets with each having at most 2/3 of the original vertices (2/3-split). Further- 
more, this separator can be found in time O(ulogn). Due to lack of space the entire algorithm 
and its analysis is presented in Appendix A. We can summarize these results as follows. 

Theorem 4. A n 1_1 / k e~ k+1 ^-separator decomposition with 2/3-splits of spanner G', for unit 
ball graphs in M k can be found in time O(nlogu). 

4.1 Proximity Problems on Spanner 

Gao and Zhang |GZ05| discuss many approximate proximity problems for the case of points in 
plane. Using a well separated pair decomposition for unit disk graph on u points in the plane, 
they show how to answer (1 +e) approximate distance queries in time 0(1 ) after O (n^/ n log ne~ 3 ) 
time preprocessing. They also estimate the diameter of the point set in time 0{n^T\ logne~ 2 ). 

We now note that other than the standard advantages of using a sparse spanner for solving 
approximate proximity problems, the separator decomposition of the spanner buys us a little 
more. We briefly describe these advantages for ball graph G in M k with f (n, e) = n ' k e — k+1 ' 2 . 

• 3/2 + e-Diameter Approximation: Due to lack of space, the algorithm is presented in 
Appendix B. The running time is of 0(uf (u, e) + f 3 (n, e)). 

• Distance Queries: It is easy to extend the results of shortest distance algorithms for 
planar graphs from Arikati et al. |ACC + 96 to spanner G'. By using the same technique, 



after 0(nf (n, e)) time and space preprocessing, we can answer an (1 + e)-stretch query in 
0(f(n, e)) and a (2 + e)-stretch query in time 0(logf(n, e)). 



5 Spanners for Disk Graphs 

In this section we extend the spanner results for arbitrary disk graphs. Each disk is associated 
with a level, a disk D p is of level I if: 2~ 1 < r p < 2~ 1 . Again we first present a high-level 
description of the algorithm before describing the actual implementation. We now stress the 
changes in construction for disk graphs from the case of unit disk graphs. Let G denote the disk 
graph and G' be the constructed spanner. 

Constructing the Forest: Let V denote the forest from the quad-dissection of the e-grid. We 
introduce a disk of level I only at depth I. That is for a node t € T of depth dt, the set C(t) 
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consists of only disks which are of at most level dt and their centers lie are contained in node t. 
The definition of an interesting node remains the same. 

Choosing the Representatives: For leaf nodes representatives are chosen in the same manner. 
For an internal node t, pick a child t' with R t / G C(t) and set Rt to R t /. 

Neighborhood of Nodes: A node t belongs to Roots if: (i) d t = and C(t) 7^ 0, i.e., square 
of e-grid with non-empty C(t), or (ii) C(t) ^ and C(P(t)) = 0. For every interesting node t, 
define Nc(t) as the set of all nodes at depth d t which are within a distance of 2~ dt from t. For 
every interesting node t G' Roots, set its far neighborhood Np(t) to 0. 

For a node t G V and integers a, |3 define its [a, (3] -shift, 
(S([a, |3],t)) with a G [-2e- 1 ,-(2e)- 1 ]U[(2e)- 1 ,2e- 1 ] and |3 G 
[— 2e _1 , — (2e) -1 ] U [(2e) _1 ,2e _1 ] as the square t' obtained by 
shifting the x- and y- coordinates of all points in t by ea2~ dt 
and e|32~ dt respectively. For every t G Roots, construct 0(e~ 2 ) 
ordered buckets, where bucket ([a, (3],t) is defined as 

bucket( [a, (3] , t) = {S ( [a, |3], t) , S( [a, (3] , P(t) ) , S ( [a, |3] , P(P(t)) ) , . . .) 

Empty squares can be ignored. All nodes in a bucket are also actual nodes in V. For t G Roots, 
do the following: scan through every bucket ( [a, (3] , t) and find the first node t' G bucket ( [a, (3] , t) 
such that there exists v G C(t) and u G C(t') with edge (u,v) in G Add t' to Np(t). 
Edges in Spanner G': Added exactly the same manner as with unit disk graphs. 

Lemma 6. The spanner G' has 0(ue~ 2 ) edges and is a subgraph of disk graph G. 

Proof. Same as Lemma El □ 
The idea behind creating the buckets is (i) to ensure that there exists lines passing through all 
the nodes in a single bucket and t, and (ii) to ensure that disks that intersect any disk centered 
in t, have their centers either close or within a node of some bucket. The first fact easily follows 
from the construction. Before proving the second fact it is convenient to define for a node t G P, 

Bucket(t) = |J bucket([a, (3],t) 

a,|3 

with a, (3 in the interval [-2e _1 , ~(2e)" 1 ] U [(2e) -1 ,2e -1 ]. 

Lemma 7. Let u be a disk center contained in node t. Then for every edge (u,v) in G, v is 
contained in Bucket(t) U Nc(t). 

Proof. If d(u,v) < 2~ dt then v is contained in Nc(t). If 2~ dt < d(u,v) < 2~ dt+1 then v is con- 
tained in U«,pS( [a, p],t). If2- dt+1 < d(u,v) < 2~ dt+2 then vis contained in U a> pS( [a, |3],P(t)). 
Similarly for increase in the distance by a factor of 2, we need to move one position in all the 
buckets of t to finish the proof. □ 















5],t) 


u 

















S([5,5],P(t)) 



Lemma 8. Let (u, v) be a edge in the disk graph G. Then there exists a path in the spanner G' 
such that dc(u,v) < (1 + c,je)d(u,v) for some constant c^- 

Proof. See Appendix C. □ 



9 



5.1 Algorithmic Version 

Constructing the Forest: We first construct the forest V exactly as described for unit disk 
graphs. At every node we also maintain the disk of the largest radius in it. The representative Rt 
of t G T' is chosen to be the disk of the largest radius, which was also the representative in one of 
the children of t. Then we visit each node t and remove the edge to its parent P(t) if the largest 
radius among disks with centers in t is less than 2~ dt . If the edge is removed then t is added to 
Roots. 

Finding the neighborhoods: The close pseudo-neighborhood of a node can be obtained as 
with unit disk graphs. For determining the far neighborhood of a node t G Roots we make the 
following distinction. 

Small Global Stretch: In this we use the observation from Gupta et al. GJS99] that maps the 
problem of reporting the intersection of a collection of disks with a query disk into halfspace range 
searching in two dimensions higher. Aggarwal and Matousek AM!)") construct a data structure 
for R k , which for any parameter n < m < nL k / 2 J and any positive constant 6, after 0(m 1+5 ) space 
and time preprocessing answers half-space queries in time 0(n/m 1// '- lc / 2 J) ) and has 0(m 1+6 /n) 
amortized update time. We construct this data structure bottom-up. Let t be a node in V with 
children's tt,i = i, ...,4 satisfying |C(tj,)| < | C (tzj.) | . We construct the data structure for t by 
inserting C(ti), C(t2) and C(t3) into C(t4). We then query the data structure with every disk in 
C(t') with t G Bucket(t') to check for intersection. 

If p{T>) is the global stretch, then each bucket([a, |3],t) is of size 0(log(p(D))). Therefore 
|Bucket(t)| = 0(log(p(£>)e~ 2 )) and Bucket(t) can be constructed within the same time bound. 
Each disk acts as a query disk log(p(2?))e~ 2 times, so total number of queries is at most O(ulog 
(p(P))e~ 2 ) times. 

To balance the total time for setting up the data structure at every node and the total query 
time, we assume the parameter m of AM95 to be n c for some c > 1 . As we work in two 
dimensions higher, each query can be answered in time 0(u 1_c / 2 ). The total time for answering 
queries is 0(u 2 ~ c / 2 log(p(D))e~ 2 ). The total time for setting up all the data structures of |A~M95 
by the procedure described above is 0(n c+6 ), where 6 is any positive constant. On eliminating 
the constant c, by balancing the exponents of u in these two running times, we get a running 
time of 0(u 4 / 3+5 log 2 ^ 3 (p(P))e -4 / 3 + ue~ 2 ) for finding the far neighborhood of all nodes. The 
space used for storing the data structures is 0(u c+5 ). 

Large Global Stretch: In this case we use adjacency lists to place the edges (u, v) satisfying 
(i) u G C(t), v ^ C(t), and (ii) r u < r v in the right bucket. Consider the line segment (edge) 
connecting u and v, and find points v' and v" on the line segments such that 2d(u,v') = d(u, v") 
and v' is contained in (J a ^ S([oc, |3], t) and v" is not contained in (J a ^ S([a, |3], t) with the usual 
ranges of a and |3. Let a', (3' be such that v' is contained in S[[oc', |3'],t). We put edge (u,v) in 
bucket([a', |3'],t). Once all the edges have been put into their respective bucket, we pick the 
shortest edge in each bucket and add it to the spanner. The entire procedure can be implemented 
in 0(|E|) given the adjacency list. 

5.2 Extension to Ball Graphs 

The extension to higher dimension can be carried out as described for unit ball graphs. For balls 
graphs in IR k , the size of N^ft) = 0(e~ k ) and Nj=(t) = 0(e~ k ). Again we can save a factor of 
by sparsifying as for unit ball graphs. In R k using results of range-searching from AM9&j we 
get a running time of 0(n 2d+6 log d p(D)e~ kd + ne~ k ), where d = 1 — 1/(|_k/2J +2). For large 
global stretch the running time is 0(|E| + ue~ k ) as we use adjacency lists. We can summarize the 
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results as 



Theorem 5. There exists an algorithm for producing a (1 + e) -spanner for ball graphs in M k . 
The spanner has 0(ne~ k+1 ) edges and for any positive constant 6 can be constructed in time 
0(mm{n 2d+6 log d p(P)e- kd ,Adj(P)} + ne- k ), where d = 1 - 1/(|k/2J +2) and Adj(D) is the 
time for constructing the adjacency list for V. 

6 Concluding Remarks 

A special case of disk graphs is when all disks have almost the same radius (between [1 — e, 1]). In 
that case our spanner construction for unit disk graphs still works by using Sharir's |Sha85j red- 
blue disk intersection algorithm instead of bichromatic closest pairs. Therefore we get a running 
time of 0(nlog 2 ne~ 1 +ne~ 2 ). Additionally all our results on spanners also extend to cases when 
intersections are between squares, regular polygons, or congruent axis-parallel rectangles, as well 
as their higher dimensional versions. 

The bottleneck in the running time of our algorithms is the time for finding the intersections 
between sets of balls. Faster algorithms for finding the intersections would improve the running 
time of the construction. 
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7 Appendix A: Separators in Spanner Graph for Unit Ball Graphs 

Let G' be the spanner for unit ball graph G constructed in Sectional We now show that the G' 
has a -y/ne~ 3 / 2 - vertex separator and can be found in O(nlogn). Let T{\) be the set of nodes in 
T' at depth I. Let S[\) be the set of all edges that were added to G' when the nodes in T{1) were 
considered. Let T>{\) be set of disk centers u with (u, v) G S{\). 

We use recursive partitions of rectangles. Let X be the list of disk centers sorted by x- 
coordinates. Let X{TZ) be vertices in rectangle TZ sorted by their x-coordinates (similarly we 
define y and y{TZ) for y-coordinates). We say a vertex crosses a line segment if any edge incident 
on it in G' crosses the line segment. Define a rectangle to be active if at least 2n/3 of the disk 
centers are inside it and not more than 0(^ l /^\e~^) vertices inside it crosses the boundary. The 
initial rectangle is the bounding box (and is active). 

A vertical double line separators of an active rectangle is a set of at most two vertical line 
segments that partition the active rectangle, such that only Ofy'ue -3 / 2 ) vertices of G' crosses 
each line segment (similarly define horizontal double line separators). Our algorithm recursively 
partitions an active rectangle alternatively with a horizontal or a vertical double line separators 
and stops when there is no more active rectangle. 



12 



In order to find the line separators we maintain the following data 
structure. We partition all the x- coordinates in X into intervals where 
every interval has at least y/ne^ 3 "/ 2 and at most consecutive 
points in X. The length of the interval bq, Xj] is the difference between I 
the coordinates of its first and last element, i.e., |x| — x,|. We maintain f 
a balanced binary search tree where leaves represent these intervals s 
and each internal node stores a splitting value between its left and 
right subtree. Additionally with each node we store the largest length 
in its subtree. Similarly we construct a tree for y. 
Constructing double line separators for level I: Let TZ be the currently active rectangle. 
Let X{TZ) = {xi , X2, . . .}. We maintain sorted doubly linked lists for both the x- and y- coordinates 
with pointers between elements representing same vertex point in two lists. Therefore in going 
to TZ from previous active rectangle TZ' , the lists X{TZ) and y{TZ) can be constructed in time 
proportional to number of points removed from TZ' to 1Z. 

Let x m be the median in X[TZ). We first find the leaf node containing the median. We then 
walk up the tree until we find the first left parent having an interval of length greater than 2~ l+2 
in its subtree. We access the rightmost leaf in the subtree representing the interval [xj,Xj/] with 
l*j ~ x y\ > 2~ l+2 . Define Li to be a vertical line segment in R to the left of xy at a distance 2~ l+1 . 
Symmetrically we construct L2 by considering only right parents. 

For updating the intervals when we go from 1Z' to 1Z we use two operations: merge and split 
and merge. When an interval becomes smaller than \/rie~ 3 ^ 2 we merge it with its left or right 
interval whichever is smaller than 2-^/ue~ 3 ^ 2 . If neither satisfies the condition we split one of them 
into two halves and then do the merge. Both of these operations and the subsequent updates to 
the tree can be performed in O(logn) time and are only done when an actual vertex is deleted. 

Li and L2 divide 1Z into at most 3 rectangles, of which at most one forms the new active 
rectangle. In this new active rectangle 1Z" we repeat the same procedure with y{1Z") to find 
horizontal line segments L3 and L4. The active rectangle thus formed is associated with the level 
following I in V . 

Lemma 9. The distance between I_i and L2 at depth I is less than 2 

Proof. Starting from the median we move I_i (or L2) only because there are -y/ue -3 / 2 disk 
centers to the left (or right) within a distance of 2~ l+2 . As there are only u disks, the result 
follows. □ 

Remark. The same distance bound holds also for the separation between L3 and L4. As a conse- 
quence of the above lemma we also get that length of Li or I_2 for I > 1 is at most 2~ lv+2 e^/ 2 \/r\. 
where l p is the level preceding I in V. We now show that Li and L2 form vertical double line 
separators for TZ. The arguments for L3 and L4 are similar. 

Lemma 10. Let TZ be an active rectangle of depth I. Li and L2 define a vertical double line 
separators for TZ. 

Sketch of Proof. We work with Li . The case for L2 is symmetric. Since TZ is an active rectangle 
there are only 0{^/rie~ 3 '/ 2 ) vertices, having edges of type (u,v), u ^ X{TZ) and v £ X{TZ) that 
cross TZ, and possibly L-|. In the remaining cases we only consider edges with both endpoints in 
TZ. The edges in Ui'>i^(t') have length at most 2~ l+1 . Any edge (u,v) G Ui/>i crossing Li 
must have both u and v be within a distance of 2~ l+1 from L-|. This implies that the number of 
edges in Ui'M^U') an d hence the vertices in Ui'>i^(t') that cross Li is 0{y/ne^^ 2 ). So if I = 
we are done. 




vertical 
«->double line«-> 

separators 
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We finish the proof by showing that not more than 0{y/rie" ') vertices in Ui'<i^ > (^') cross 
Li- Fix any level I' < I. From Lemma 03 and the remark above we know that the length of I_i 
is less than 2~ l p +2 - v /ue 3 / 2 . In 1Z consider the rectangular region between L] and a line segment 
drawn at distance 2 -1 +1 from L-|. The area of this region is 2~ ^ +2 2 — * +1 e 3 / 2 -^/n and any edge of 
£ [V] that crosses L-j will have one endpoint in it. The number of nodes of T{V) in this region is 
2 l — lp+5 ^/rie^ ^ 2 . For each node of T(l') at most e — 1 vertices contained in it could be present in 
T>{\'). Thus the number of vertices in D(l') in TV that cross Li is at most 2 l ~ lp+3 A/n~e~ 3 / 2 . This 
forms a decreasing geometric series in I'. For levels I" which are more than O(logn) below I, 
only constant number of disks of T>{\") crosses I_i . Our method of picking representatives ensures 
that they are the same over all such I". Summing over all I' < I we get that the total number of 
vertices in Ui'-ci^U') crossing I_i is 

Y_ 2 l '- l P +3 v^e- 3/2 = 0(v^e- 3/2 ) 

0<l'<l p 

□ 

Final Shape of Separator: We stop the algorithm if no rectangles are active. The rectangle 
with the largest number of disk centers form one component of the separator. The vertices of it 
which are connected to the outside are the vertex separators. 

The proof of termination is omitted in this extended abstract. The algorithm also extends 
to higher k and we get a n /^e - k+ ' 2 -separator, The rectangle gets replaced by k-dimensional 
box and line segments are (k — 1 )-dimensional hyperplanes. We recursively partition along every 
dimension. 

Appendix B: Diameter Approximation 

We use a rooted binary tree Jq to represent a separator decomposition of G = (V, E). For a set V 
of vertices in G, we use N(V') to denote the neighborhood of V. Each node t G Jq is labeled by 
two subsets of vertices V(t) C V and S(t) C V(t). Let G(t) = (V(t),E(t)) denote the subgraph 
induced by V(t). Then S(t) is the separator in G(t). The root r G Jq has V(r) = V and S(r) is a 
separator in G . For any t G Jq, the labels of its children to, ti are defined as follows: Let Vi C V(t) 
and V 2 C V(t) be the components separated by S(t) in G(t). Then V(t ) = V, U (S(t) n N(Vi)), 

v(ti)=v 2 u(S(t)nN(v 2 )). 

We describe our 3/2-diameter approximation algorithm for any graph G = (V, E) from the 
family of weighted directed graphs having a M ^-separator decomposition, i.e., every subgraph 
of size M has an M^-separator. The algorithm is similar to the algorithm for 3/2-diameter 
approximation of Aingworth et al. ACIM99] for general graphs. On graphs with m edges their 
algorithm runs in time 0(m.\/n logn + n 2 logn). 

For e G E, let w(e) denote the weight of edge e. We associate boundary vertices, B(t) with 
each node t G Tq. The boundary of the root r G is B(r) = 0. The boundary of every other 
node t is defined as B(t) = S(P(t)) U B(P(t)) n V(t), where P(t) is the parent of t in Jq. We 
discuss the procedure for single node t G Jq. 

Performing SSSP: Do Single Source Shortest Path (SSSP) from S(t) on the graph H(t) = 
(V(t),E(t)UB(t)xB(t)), with weights w H defined as follows: For a edge e = [v h v 2 ) G B(t)xB(t) 
define Wu t )(e) to be d G ( t )(vi , v 2 ). One can inductively see that these weights are available. For 
t = r it is true. If vi or v 2 is in B(P(t)), inductively we know that weights are available. Otherwise 
both vi and v 2 are in S(P(t)) and we get the weights after doing SSSP on H(P(t)) (see Lemma 
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ITTjl . For edge e G E(t), w H(t) (e) = w(e). 

Estimate for a node: Let maxi be the longest distance obtained in H(t) after performing the 
SSSP from all vertices in S(t). For u G V(to) U V(t)), let clo(u) denote the closest vertex in S(t). 
Let z be a vertex with clo(z) = max{clo(u)}. For z G V(tj), do SSSP from z on H(t) to find a 
farthest node w G V(t-j_j), j = 0, 1. Let max2 be the distance between z and w. Compute the 
estimated diameter for node max(t) as max{maxi , max2}. 

Output of Algorithm: Output the estimated diameter Dia as max{max(t) : t G Tq}. 
Lemma 11. For any t G Tq and u,v G V(t), dH{t)(u-,v) = dc(u,v). 

Proof. Proof by induction. If t = r, the claim is true. By the inductive hypothesis dn{p(t)) ( u > v ) = 
dc(u,v). Let P = {u,pi,p2, . . . ,v} denote the shortest path from u to v in H(P(t)). If P C V(t) 
then the shortest path is also present in H(t). Otherwise let i be smallest index such that 
Pi G S(P(t)) and let j be the largest index such that p, G S(P(t)). Both pt and p, are in B(t) 
and have an edge connecting them with weight w H ( t )(p-t,pj) = dctp-bPj)- Thus distances are 
preserved in H(t). □ 

Lemma 12. Let G be a graph with n vertices and an M^- separator decomposition, < u < 1, 
with ct-splits, then number of edges in G is O(n). 

Proof. The number of edges E(n) in G satisfies the recurrence: E(n) < max a {E(an + Cn^) + 
E((1 — a)n + Cu^)}, where C is a given constant from the size of the separator. We can show 
by induction that E(n) < cn — dn^ for constants uo,c and d chosen as follows : (i) choose 
n-o such that u — (2Cn^ i )/(a ,x + (1 — a)^— 1) > 1 for all n > no, (ii) choose c = ), where 
ui = max{no/a, no/(1 — a)}, and (iii) choose d such that 2C/(a (X + (1 — a)^— 1) < ~. □ 

Theorem 6. Let G be a positive weighted graph with an M^- separator decomposition. The output 
Dia provides a 3/2 approximation for the diameter in G and can be computed in time 0(n 1+ ^ + 

Proof. Let a and b be the two vertices whose shortest path defines the diameter. Let dc(ci,b) = 
A. Let t be the node in Tq with a G V(to) and b G V(t-j), where to,ti are the children of t in 
the separator decomposition tree. If either a G S(t) or b G S(t) our estimate for diameter equals 
A (from Lemma ITT|l . We now consider two cases: 

Case 1: Vu G V(t ) UV(ti), 3v G S(t) such that d G (u,v) < A/3: Let w be the node in S(t), s.t. 
dG(a,w) < A/3. Then definitely de(w, b) > 2A/3 . This implies that we get 3/2 approximation 
for the diameter as we do SSSP from w. 

Case 2: Otherwise: At least one node exists in V(to) U V(t-|) to which the distance from all 
the nodes in S(t) is at least A/3. This implies that z G V(to) U V(ti) is also at least A/3 from 
all the nodes in S(t). Assume z G V(t-j). If the farthest node from z is at a distance 2A/3, 
this distance is our estimate and we are done. Otherwise, let pi G S(t) be a node through 
which the shortest path from z to a passes. dc(z, a) = dc(z,pi) + dc(pi, a) < A. Since 
dc(z,a) < 2A/3 and dc(z,pi) > A/3. This implies dc(pi,a) < A/3. Since dc(a,b) = A, 
implying dG(b,pi) > 2A/3. Again we get a 3/2 approximation when we do SSSP from p-|. 

To analyze the running time, note that for all t G T^, |B(t)| = O(n^), therefore at every t we 
introduce at most 0(n 2(X ) edges. The total running time is dominated by the cost of doing SSSP. 
Using Dijkstra for SSSP and Lemma IT21 gives the required bounds. □ 
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Appendix C: Proof of Lemma [S] from Section [5] 

Proof. Assume w.l.o.g. that r u < r v . Let a be the deepest node in Roots containing point u. 
Let t be the node containing v with d Q = dt- We divide the analysis into parts based on the 
position of t. For the case when t S Nc(a) the arguments are similar to those in Lemma 0] and 
one can show d G /(u,v) < (1 + c u e)d(u,v). 

If t ^ N G (a). Consider bucket([a', (3'], a) containing v. From Lemma0we know such a', (3' 
exists (as we are in the case where v is not contained in Nc(t)). Let b G bucket([a', a) be a 
node containing v. Let f € bucket([a', |3'], a) n Np(a). Let (u',v') with u' G C(a) andv'e C(f) 
be the edge added to the spanner. The edge (u, u') also exists in G because the the minimum 
radius disk in C(a) is at least e _1 times a's length. In G' we know from the previous case that 
d G /(u,u') < (1 + Cu.e)d(u, u'). 

Let ai and az be the ancestors of a in V at depth 
db and df respectively. Let f 1 be the ancestor of f in V 
at depth db- In the following description the distance 
between two squares is the distance between their cen- 
ters and distance between a point and a square is the 
distance of the point from the center of the square. 

We first claim that d G /(u,v) < (1 + Ci e)d(a-j , b) 
for some constant c-\. Since f is obtained by shift- 
ing az by ecx'2~ da 2 and e|3'2~ da 2 , we get d(u, az) < 
c s ed(a2,f), d(f,v') < c s ed(ci2,f), and d(u.u') < 
c s ed(u, v'). Applying the triangle inequality we ob- 
tain 

d(u,v') < d(u,a2) + d(a2,f) + d(f,v') < (1 +2c s e)d(a 2 ,f) 

In G' we have 

d G /(u,v') = d G '(u,u') + d(u',v') < (1 + c u e)d(u,u') + d(u',u) + d(u,v') < (1 + c 2 e)d(u,v') 

Substituting for d(u,v') implies that d G /(u,v') < (1 +Cie)d(aa,f) for some constant c-\. Note 
that d(a 2 ,f) = d(ai,fi). To complete the claim we use induction over distances. Since the edge 
(v',v) also exists in the graph by the inductive hypothesis we get d G /(v',v) < (1 + Ci e)d(f-| , b). 
By construction d(ai,b) = d(ai,fi) + d(fi,b), therefore we get d G '(u,v) < (1 + Cie)d(ai,b). 

To complete the proof we repeat similar arguments to show that d(a-|,b) < (1 + c u e)d(u,v). 
Substituting for d(ai,b) in terms of d(u,v) completes the proof. □ 
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